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ABSTRACT 

We present results based on Spitzer Space Telescope mid-infrared (3.6-30 /im) observations of the 
nearby IIP supernova 2005af. We report the first ever detection of the SiO molecule in a Type IIP 
supernova. Together with the detection of the CO fundamental, this is an exciting finding as it may 
signal the onset of dust condensation in the cjecta. From a wealth of fine-structure lines we provide 
abundance estimates for stable Ni, Ar, and Ne which, via spectral synthesis, may be used to constrain 
nucleosynthesis models. 

Subject headings: supernovae: general — supernovae: individual(SN 2005af, SN 1987A, SN 2004dj) 



1. INTRODUCTION 

Mid-infrared (mid-IR) studies of supernovae (SNe) 
have been hampered by a combination of several factors: 
their intrinsic faintness at epochs of interest, the high 
background at mid-IR wavelengths, and strong terres- 
trial atmospheric absorption. With the sol e exception of 
the extremely nearby (50 k pc) SN 1987A ijRoche et al. I 
119931: IWooden et al. 1 11993). this wavelength region has 
remained out of reach for SN studies. Yet, access to 
the mid-IR region is important for at least two reasons. 
First, a large number of fine-structure lines that are sen- 
sitive to the details of the explosion physics lie in the 
5-40 /jm region. These provide robust measurements of 
element abundances permitting direct comparisons with 
explosion models. Second, core-collapse SNe (CCSNe) 
have been pro posed to be major sourc es of dust at high 
redshifts (e.s.. iTodini fc Ferrarall200l . but even in the 
local Universe, this hypothesis rests on little direct ob- 
servational evidence. By monitoring the mid-IR evolu- 
tion, where warm dust emits most strongly, we have a 
powerful means by which we can ascertain the epoch of 
dust condensation, and estimate the amount of fresh dust 
produced in the ejecta. Alternatively, via mid-IR mea- 
surements, we can study the development of an IR echo 
arising in the pre-existing dusty circumstellar medium, 
thereby constraining the pre-supernova evolution of the 
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progenitor (e.g.,Bj et~aT~ll2006lL 

The improvements in terms of sensitivity and spa- 
tial resolution afforded by the Spitzer Space Telescope 
(Spitzer) finally allow us to probe SN behaviour using 
mid-IR diagnostics. To this end, we are pursuing a multi- 
epoch spectroscopic and photometric campaign of a sig- 
nificant sample of SNe of all types, within the frame- 
work of the Mid-IR Supernova Consortium (MISC). For- 
tuitously, since the launch of Spitzer, at least half a dozen 
CCSNe have exploded within ~10Mpc. Recently, we 
presented a first mid- IR study of th e nearest (3.13 Mpc) 
of these, SN 2004dj ijKotak et al. 1 12005). Here, we de- 
scribe results for the Type IIP SN 2005af that highlight 
the diversity of this subgroup. 

SN 2005af was discovered in the nearby edge-on spi- 
ral galax y, NGC 4945, on 2005 Feb. 8.22 (UT) at V= 
+12.8 bv l.Iacaues fc Pimentel I il2005j). On the basis of a 
spectr um obtained four days later. iFilippenko fc Folevl 
( 2005) classified it as a Type IIP (plateau) supernova, at 
an epoch of about one month post-explosion. In what 
follows, we assume that the supernova exploded on 2005 
Jan. 9, but note that the explosion epoch is probably 
uncertain by up to a few weeks. NGC 4945 is associated 
with the Centaurus group of gal axies, and estimates of 
its di stance have r anged from 3 .9 l|de Vaucouleurs et al. I 
11981(1 to 8.1 Mpc ijBaan 111985]) . In recent years, these 
values have tended to converge ne ar the lower end of this 
range. Following the reasoning in lBergman et aTl <|1992l 
and references therein), we adopt a distance of 3.9 Mpc. 

2. OBSERVATIONS 

We first observed SN 2005af with Spitzer during Cy- 
cle 1 using the Infrared Spectrograph (IRS) and IRS 
Peak-up Imaging (PUI) under Director's Discretionary 
Time Program 237. We subsequently observed it dur- 
ing Cycle 2 as part of our GO Program 20256, using the 
Infrared Array Camera (IRAC), IRS, and IRS-PUI. 

2.1. Photometry 

3.6-8.0 /urn imaging with IRAC was obtained on 2005 
Jul. 22 (194 d) and 2006 Feb. 12 (399 d). A frame time 
of 30 s and 5 (medium-sized) dither positions in mapping 
mode for each pair of channels were used. Our mea- 
surements from aperture photometry, carried out on the 
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Fig. 1.— Mid-IR spectra of SN 2005af for days 67 (top curve) 
and 214 (bottom curve). Plausible line identifications are marked. 
The IRAC measurements are shown as triangles (194 d) and (red) 
diamonds (399 d), while the black squares show the PUI measure- 
ments for 67, 211, 433 d from top to bottom, respectively. Note the 
sharp drop in flux between the two IRAC epochs at 4.5 fim, and 
the relative increase in prominence of the flux at ~8 fim. The IRS 
flux levels are consistent with the photometry. The insets show the 
evolution of line profiles described in J3] 

post-Basic Calibrated Data, vS14.0.0 image mosaics us- 
ing GAIA, are shown in Table □ At 67 d, SN 2005af is 
very bright in all four IRAC channels, The flux is partic- 
ularly bright in channel 2 (4.5 fim). Nevertheless, exami- 
nation of individual BCD frames showed that saturation 
had not occurred i.e. the measurements are robust. 

Our first epoch of PUI-16/im imaging was obtained 
in Peak- Up mode only, i.e., before dithering capabilities 
were offered with IRS. The two subsequent epochs of PUI 
observations were carried out using 31.46 s ramps and 20 
small-scale dither positions. The SN was clearly detected 
at all three epochs (Table QJ. 

2.2. Spectroscopy 

Spectra were obtained on 2005 Mar. 17 (67 d) and 
2005 Aug. 11 (214d) with the Short-Low (SL, A = 5.2 - 
14.5 fim; R = 64-128) module of the IRS in staring mode. 
A Long-Low (LL, A = 14-38^m; R = 64-128) spectrum 
was also acquired at 67 d. For the SL mode, the total 
time-on-source at the two epochs was 1219 s and 1828.5 s 
respectively, and 629 s for the LL mode. 

The first and second epoch data were preprocessed 
using versions S12.3.0 and S12.4.0, respectively, of the 
Spitzer data processing pipeline. We began the reduction 
of the BCD data by differencing the two nod positions to 
remove the background emission. All subsequent steps, 
i.e., straightening, extraction, wavelength and flux cali- 
bration were carried out using SPICE vl.l-betal5 with 
standard settings The spectra are shown in Fig. 2] We 
also extracted the spectra using post-BCD S13.2.0 prod- 
ucts and found that they were virtually identical. 

3. ANALYSIS 

The high IRAC flux at 4.5 fim for the first epoch com- 
pared to other channels is almost certainly due to CO 
fundamental band emission. This is corroborated by the 
steep rise towards 5 fim in both spectra, which we identify 
with the red wing of the CO fundamental (see Fig. 0. 
Between 194 d and 399 d, the sharpest decline in IRAC 
flux is seen at 4.5 fim. These observations suggest that 



CO formed in the ejecta less than 2 months after ex- 
plosion, persisting to at least 7 months, but had faded 
substantially by about 1 year. The relative increase in 
the IRAC channel 4 (8 fim) flux by 399 d is probably due 
to the emergence of spectral features in the 6.43-9.40 fim 
region spanned by this channel. This is discussed below. 

By day 214, a pedestal-like feature appears spanning 
~7.6-9 fim, which we identify with molecular emission 
due to the SiO fundamental band (Figs. ^13- This is 
the first time that SiO has been detected in the ejecta of a 
SN IIP the most common subtype of all CCSNe at low 
rcdshifts. The only other detection of Si O in SN ejecta 
was i n the peculiar Type II SN 1987A ijRoche et al. I 
ll991tlw"ooden et al. 11993^ . That molecules are detected 
at all, provides a diagnostic of the prevalent conditions in 
the ejecta. As both CO and SiO are powerful coolants, 
they are regar ded as precursors to d ust condensation in 
the ejecta (e.g. iNozawa et al~1l2003|h 

In Fig. [5] we compare a n LTE model SiO spe ctrum 
for SN 1987A at 260 days (|Liu fe Dakarno l ll994) with 
the day 214 SL spectrum of SN 2005af. The model, at- 
tenuated by a factor of 2.5, produces a convincing re- 
production of the observed behaviour. Scaling from the 
model, the ma tch to our spectrum sugg ests a SiO mass of 
~2xlO- 4 M Q . iLiu fc Dalgarnol fj9Mll996|) argue that 
such a mass can only form if there is no mixing of helium 
into the oxygen-silicon region of the ejecta. 

Apart from molecular emission, several conspicuous 
lines are visible in both spectra (Fig. w hich were iden - 
tified by analogy with SN 1987A ijWooden et al. l fl993'l. 
At both epochs, the spectrum is dominated by forbidden 
lines of Ni, Co, and Ar, of which the [Nill] 6.64 fim, [Aril] 
6.99 fim, and [Coll] 10.52 fim are the most prominent. 

Between 67 d and 214 d the lines show appreciable in- 
creases in intensity (see Table [2 Fig. QJ, In addition, in 
spite of the low spectral resolution (R=64-128), we ob- 
serve in a number of lines (e.g., [Nill] A6.64/im, [Coll] A 
10.52 fim) an extended red wing at 67 d which fades by 
214 d. The broad blend with an extended red wing near 
11.3 fim at 67 d is probably due to a combination of [Nil] 
and HI (7-9), but by 214 d, it is symmetric. Extended 
red wings were also detected i n several lines in SN 1987A 
Ce.g.. lWitterborn et al. 11989]) accompanied by small red- 
shifts of the line centroid. A plausible explanation is pro- 
vided by electron scattering in the exp anding hydrogen 
envelope ( Fransson fe Chevalier I [1989(1 . This may also 
apply to SN 2005af, although we do not detect any ob- 
vious redshifting of line centroids. 

We now estimate the abundances from the stronger 
lines at the later epoch when the ejecta has a lower 
optical depth. Following iMcCrav I Jl9901. we estimate 
the electron density of the ejecta as n e ~ 10 8 i~ r 35 

cm~ 3 w 7 x 10 8 cm~ 3 at day 214. We derived critical 
densities for Ni + and Ar + u sing collision strengths from 
IBautista fc PradhaiTI <J1996|) and IQsterbrock fc Ferland I 
l|2006[K respectively, and found these to be 1-3 orders 
of magnitude smaller than n e . On this basis, we deem 
an LTE treatment to be valid. Using the line intensi- 
ties in Table [3 we derive masses using the upper ex- 
cited state populations of the ions and find 3.7xl0~ 3 Mq, 
1x1O~ 3 M , and 2.2xlO~ 3 M Q for stable (Ni°,Ni+), 
Ar + , and Ne + respectively, for an electron temperature 
of 4000 K. Note that the result is insensitive to temper- 
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Fig. 2.— The 260 d LTE model from lLiu fc Dalgarnol JTflflHV 

attenuated by a factor of 2.5, superimposed on the continuum- 
subtracted second-epoch spectrum of SN 2005af. From this, we 
deduce a SiO mass of ~ 2 x 1CT 4 M (©. 

ature: a ±1000 K change in temperature, changes the 
mass estimates by <15%. The partit ion functions are 
from lHalenka et al. 1 1)2001() for Ni, and llrwin 1 l)1981|) for 
Ar and Ne. Assuming a uniform distribution of stable Ni, 
we find that the Sobolev optical depths of all detected 
Ni lines are less than 0.06. Strong transitions of [Nilll] 
11.002 fim and [NilV] 8.405 /im lie within the wavelength 
region covered, but these lines are not detected at ei- 
ther epoch. We conclude that the Ni is predominantly 
in the neutral and singly-ionised state. Thus, assuming 
insignificant clumping, the value we derive above is prob- 
ably close to the total stable Ni mass. 

The predicted mass of stable Ni is sensitive to the 
progenitor mass, although this relation is not mono- 
tonic, and depends sensitively on the explosion mecha- 
nism, multi-dimensional effects and, in spherical mod- 
els, on the mass cut, which is a free parameter. 
iThielemann et al~l l)1996f ) predict stable Ni masses in 
the range 0.01-0.014 Mq for progenitors in the range 
13-20 Mq, and a much lower value of 0.002 Mm for a 
25 Mq star. Other authors (e.g.lChieffi fc Limongi 120041: 
iWooslev fc Weaver! 119951 iNomoto et al. 1 12006ft derive 
comparable values, though detailed intercomparisons are 
difficult due to different parametrizations of physical pro- 
cesses; 

iChiefE fc Limongil 1)20041) provide explosive yields as a 
function of metallicity based on parameters which have 
been calibrated to fit observed properties of massive 
stars. The value for the stable Ni mass that we derive 
is closest to their progenitor models of 13-15 M Q with 
metallicities of a third to a twentieth of the solar value. 
Although we can definitively rule out higher mass pro- 
genitors at even lower metallicities, we current cannot 
exclude a 35 Mq progenitor at solar metallicity. Further, 
planned, late-time spectroscopy should provide tighter 
constraints. 

Given the high energy required to ionise neutral Ne 
(21.6eV) and Ar (15.8eV) we conclude that, given the 
low ionisation conditions indicated by Ni, the bulk of 
the Ne and Ar lies in the neutral state. However, no 
strong transitions of [Nel] or [Arl] lie in our wavelength 
range. Consequently, the masses derived here for the first 
ionisation state represent lower limits only. 

F or SN 1987A at 260 days, using the [Coll] 10.52 ^m 
line. lRoche et al.1 1)1993(1 derive a value for M Co + of 4.6 x 
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Fig. 3. — Comparison of continuum-subtracted spectra of CC- 
SNe at roughly coeval epochs. The dotted lines mark the zero- 
level. Note the striking difference in the strength of the [Aril] 
line at ~200 d between the SNe. The SN 1987A spectra are from 
IWooden et al. I J1993I) ; the 260 d spectrum has been divided by 
a factor of 5. Isolated features that are common to both epochs 
are labelled in the top panel only; features that are evident in the 
later epoch only are indicated in the lower panel. The SN 2004dj 
spectrum is part of our ongoing monitoring programme; a fuller 
analysis of these data will be presented elsewhere. 

10~ 3 M Q . Using our day 214 intensity for this line in 
SN 2005af we obtain M Co + = 2.4 x 1O" 3 M . Allowing 
for the different epochs for the two SNe, this implies that 
a 56 Ni mass of 0.027 Mq was ejected by SN 2005af. This 
estimate f alls in the exp ected range for normal Type II 
SNe fe.g.. lHamuv1l2003l) . 

Although we have a clear detection of the SN 2005af 
spectrum between 16-30 /im (Fig. ^) - the first, for any 
SN at this early epoch - we do not see any significant 
spectral features. This is not surprising, given its youth. 
In SN 1987A, lines of [Fell] in this spec tral region were 
only detected after ~400d l)Haas 111990)) . 

3.1. Comparison with other Type II SNe 

In Fig. |21 we compare roughly coeval mid-IR spectra 
of SN 2005af with that of SN 1987A (Ilpec) at ~60d, 
and with those of SN 1987A and SN 2004dj (IIP) at 
~250d. The 60 d mid-IR spectrum of SN 1987A is dom- 
inated by HI lines ( see Fig. |3] or iRoche et al. 1 119931 
IWooden et al. lll993)) . In contrast, strong lines of [Nill] 
and [Aril] are present in SN 2005af. This suggests ei- 
ther that our first-epoch SN 2005af observations were 
taken on the radioactive tail of the lightcurve and that 
the length of the plateau phase must have been <60d, 
or that the estimated explosion date, based on the clas- 
sification spectrum is too late by >1 month. 

At the later epoch al l three SNe exhibit pro minent lines 
of nickel and cobalt. IWooden et al. I l)1993l) emphasize 
that the ratio of [Nil] 7.5 /im to [Nill] 6.64 /im lines pro- 
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vides a sensitive measure of the ionization fraction of Ni 
(XNi = "Ni+/ n Ni)- Using their Fig. 9, and assuming 
LTE, we find XNi ~ 0.35,0.55, and 0.85 for SN 2004dj, 
SN 2005af, and SN 1987A, respectively, at -250 d. 

SN 2005af shows much stronger [Aril] than SN 1987A. 
Moreover, SN 2004dj is found to be a factor of 20 weaker 
than SN 2005af in this line. Even if we normalise to 
the strengths of the adjacent [NiII]6.63 /mi, [Nil]7.50/tm 
lines, the [Aril] line is still over x5 stronger in SN 2005af. 
This is intriguing given that SNe 2004dj and 2005af are 
both of Type IIP. One possibility, indirectly implied by 
XNi above, is that the conditions required to ionize Ar 
are not attained in SN 2004dj, while Ni + , with its lower 
ionization potential (7.6eV), is easier to produce. 

Another hint may come from the spherically symmetric 
iThielemann et al~1 l)199fil) models, where the bulk of the 
stable Ar (mostly 36 Ar) is formed in a shell just above the 
stable Ni region. Thus, one might expect that a higher 
mass cut (below which the core collapses to a neutron 



star or black hole) within the Ni zone would produce a 
lower absolute Ni mass, but a higher Ar/Ni mass ratio. 
However, this seems to be contrary to what we observe, 
i.e., SN 2004dj has an apparently smaller Ni mass but an 
even smaller Ar/Ni mass ratio than SN 2005af. We sus- 
pect that to account for this behaviour, a certain degree 
of asymmetry or dumpiness, as well as mixing in the 
eject a has to be invoked. Interestingly, iPerevra et al. I 
( 2006) infer an overall asphericity of 20% for SN 2005af 
from optical polarimetry. 

This work is based on observations made with the Spitzer 
Space Telescope, which is operated by the Jet Propul- 
sion Laboratory, California Institute of Technology un- 
der NASA contract 1407. Financial support for the re- 
search is provided by NASA/ Spitzer grant GO-20256. 
Support for this work was provided by NASA through 
an award issued by JPL/Caltcch. R.K. and S.M. ac- 
knowledge support from an ESO fellowship and from the 
EURYI Awards scheme respectively. 
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TABLE 1 
Photometry of SN 2005af 



Flux (mjy) 





Epoch 


texp 




IRAC 




PUI 


Date 


(d) 


(b) 


3.6 /tm 


4.5/tm 


5.8 (im 


8.0 /im 


16 /im 


2005 Mar. 17 


67 


1219 










5.0±0.3 


2005 Jul. 22 


194 


524 


3.96±0.02 


17.54±0.04 


9.31±0.05 


7.43±0.06 


[3.10] 


2005 Aug. 08 


211 


629 










2.83±0.03 


2006 Feb. 12 


399 


524 


0.93±0.02 


3.43±0.02 


2.24±0.02 


4.84±0.07 


[2.63] 


2006 Mar. 18 


433 


629 










2.19±0.03 



Note. — Photometry was carried out using a 5'.'5 radius aperture for all channels. This aperture was chosen to enclose the full 
diffraction pattern to the first dark ring at the PUI red limit of 18.7 /im. The aperture was positioned using WCS coordinates. Aperture 
corrections as appropriate for the outer diffraction rings were applied to IRAC channels 1-4 using Table 5.7 of the IRAC data handbook. 
No attempt was made to apply a colour correction to the measured PUI fluxes, since the choice of temperature is not obvious. However, 
we note that the correction would be <2% for hot dust (500-1000 K), and that such a correction would lie well within our errors. We 
estimated the systematic error due to variations in the background by repeating our measurements with differing annuli. The errors above 
are a combination of systematic and statistical errors. The figures in square brackets under the PUI column are indications of the expected 
fluxes coeval with the IRAC values assuming a Q-band decline rate similar to that of SN 1987A. t eX p is the on-source integration time. 



TABLE 2 
Line intensities 

67 d 214 d 

A rcst Ion Inst.Res. I FWHM I FWHM 
(/im) (kms -1 ) (kms - 1 ) (kms^ 1 ) 



6.636 


[Nill] 


2700 


3.18 


4000 


9.42 


3600 


6.985 


[Aril] 


2600 


1.61 


3800 


6.16 


3000 


7.507 


[Nil] 


2400 






4.48 


4000 


10.521+ 


[Coll] 


3400 


0.78 


4500 


1.99 


3000 


11.308* 


[Nil] 


3200 


0.25 


2400 


0.78 


2900 


12.813* 


[Nell] 


2800 






0.56 


3300 



Note. — The intensities were measured after subtracting a continuum fitted by a low-order polynomial or a spline function; the FWHM 
were derived from Gaussian fits. Column 3 lists the instrumental resolution; I, the intensity, is in units of 10 — 14 erg cm -2 s — 1 . * As this 
line is symmetric, we assumed that there was no blend with [Nill] 10.68 /jm. * At 67 d this line has a probable contribution from HI(7-9) 
but by day 214, it is symmetric; above , we assume that it arises due to [Nil] only. Interestingly, the Hl(7-9) feature was not detected in 
SN 1987A at 60 d IWooden et al. H1993fl . * The 12.8 /im feature at 214 d is likely to be a blend of [Nell] A12.813/mi and [Nill] A 12.729 /im. 
Given the strength of the [Nill] 6.636 /im line, we estimate that the 12.8 /tm blend comprises roughly equal contributions from the two 
elements. At 67d, however, given the strength of the [Nill] 6.636 (im line, the [Nill] 12.729 /im line would have been barely detectable. 



